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不 同月 球 信 力 约束 下 的 地 月 Halo 轨道 转移 轨道 设计 
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摘 要 : 针对 地 月 系 L2 点 不 同 任务 需求 下 的 低 耗 能 转移 轨道 设计 问题 ,基于 不 变 流 形 理论 与 混合 优化 技术 ， 
深入 研究 了 不 同月 球 借 力 约束 与 不 同 幅 值 Halo 轨道 的 入 轨 点 ( 简称 HOI 点 ) 对 转移 轨道 飞行 时 间 与 燃料 消耗 的 影 
响 , 给 出 了 HOI 点 选择 策略 。 首 先 结合 任务 要 求 并 考虑 月 球 引 力 影 响 ,在 月 球 借 力 点 施加 不 同 约束 条 件 , 通过 微分 
修正 算法 调整 Halo 轨道 的 稳定 流 形 ,设计 月 球 到 Halo 轨道 的 转移 轨道 。 采 用 遗传 算法 与 微分 修正 算法 相 结合 的 
混合 优化 策略 ,在 同时 考虑 地 球 停泊 轨道 高 度 、 倾 和 角 、 升 交点 赤 经 与 航 迹 角 等 多 约束 条 件 下 ,对 燃料 最 优 的 地 月 转 
移 轨 道 进行 研究 。 最 后 ,分 析 月 球 借 力 高 度 、 借 力 方位 角 和 不 同 HOI 点 对 平 动 点 转移 轨道 飞行 时 间 与 燃 耗 变 化 量 
的 影响 ,对 于 考虑 月 球 借 力 的 地 月 平 动 点 转移 轨道 设计 和 与 应 用 具有 重要 的 参考 价值 。 
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A Design Method for Earth-Moon Halo Orbit Transfer Trajectory Under 
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Abstract: To design low-energy Earth-Moon L2 Halo orbit transfer trajectory for different mission requirements, a deep 
research on the influence of different constraints to the moon gravity-assisted maneuvers and Halo orbit insertion ( HOI) 
points with different amplitudes on time of flight and fuel consumption is made, and a selection strategy for HOI points is 
proposed in this paper. Firstly, combined with mission requirements and effects of the moon gravity, different constraints 
are applied to the moon gravity-assisted maneuvers. Stable manifold of Halo orbit is adjusted through differential correction 
algorithm and the transfer trajectory from moon to Halo orbit is designed. Then, in consideration of orbit altitude , orbit 
inclination, right ascension of ascending node and track angle at the same time, the hybrid optimization strategy combined 
genetic algorithm with differential correction algorithm is adopted to study Earth-Moon Halo orbit transfer trajectories with 
optimal fuel consumption. At last, the influence of moon gravity-assisted altitude and azimuth as well as different HOI 
points on flight time and fuel consumption of transfer trajectories is respectively analyzed, which has reference value for 
design and application of Earth-Moon Halo orbit transfer trajectories. 
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( 简称 EMI2 点 ) 附近 Halo 轨道 的 应 用 价值 逐日 突 

显 。 航 行 于 EMI2 点 Halo 轨道 上 的 探测 器 能 够 观测 

随 着 平 动 点 轨道 任务 的 相继 提出 ,以 及 月 球 背 面 ”月 球 背 面 , 从 而 为 月 球 背 面 的 着 陆 任务 提供 导航 方 
在 未 来 深 空 探测 中 存在 的 重要 作用 。 地 月 系统 L2 点 。 案 。 同 时 ,围绕 EML2 点 能 够 设计 出 到 达 月 球 、 火 星 
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和 其 它 行星 的 低 耗 能 转移 轨道 ,可 作为 探测 深 空 环境 
的 最 佳 落脚 点 "-”。 因 此 ,在 平 动 点 轨道 任务 中 , 探 
测 器 采取 何 种 飞行 方式 ,消耗 多 少 燃 料 示 能 设计 出 合 
适 的 转移 轨道 一 直 是 空间 领域 研究 的 热点 问题 。 

针对 平 动 点 转移 轨道 设计 问题 ,由 于 EMI2 点 
附近 Halo 轨道 的 稳定 流 形 无 法 与 地 球 停泊 轨道 相 
交 , 基 于 不 变 流 形 设计 雾 消耗 转移 轨道 的 方案 并 不 存 
在 。Zazzera 等 上 采用 遗传 算法 与 序列 二 次 规划 的 混 
合 优化 方法 ,设计 了 多 条 不 考虑 月 球 借 力 的 EML2 点 
间接 转移 轨道 。Parker 8 基于 弱 稳 定 边界 转移 策略 ， 
通过 连接 日 地 系统 与 地 月 系统 的 不 变 流 形 , 研 究 了 低 
地 球 轨道 飞 抵 EMI2 点 Halo 轨道 的 转移 轨道 特性 ， 
但 弱 稳 定 转移 存在 着 飞行 时 间 长 等 问题 。Gordon 
则 利用 稳定 流 形 理论 与 月 球 借 力 技术 ,构造 了 不 同 幅 
值 Halo 轨道 的 两 脉冲 转移 轨道 ,虽然 降低 了 机 动 速 
度 增 量 , 但 两 脉冲 条 件 下 借 力 点 方位 不 可 调整 ,存在 
着 局 限 性 。 同 时 ,飞行 耗 时 较 长 且 大 幅 值 条 件 下 燃料 
消耗 过 多 。 文献 6 - 8] 在 Gordon © 的 研究 基础 上 ， 
给 出 了 不 变 流 形 与 月 球 借 力 飞行 技术 相 结 合 的 三 脉 
冲 转 移 轨 道 方 案 。 但 针对 不 同 幅 值 Halo 轨道 ,并 未 
分 析 不 同月 球 借 力 约束 对 Halo 轨道 的 入 轨 点 ( HOT) 
选择 策略 的 影响 ,以 及 借 力 机 动 点 三 轴 速 度 增 量 大 小 
等 特征 ,限制 了 在 工程 中 的 应 用 价值 。 

表 1 与 图 1 描述 了 在 以 往 研 究 中 ,针对 EML2 点 
Halo 轨道 的 转移 轨道 不 同 设计 方法 及 参数 之 间 的 比 
较 。 其 中 ,CR3BP 为 圆 型 限制 性 三 体 模 型 , BCRFBP 
为 双 圆 限制 性 四 体 模型 ,JPL 为 星 历 模型 ,IT 为 间接 
转移 , WSB 为 弱 稳 定 转 移 ,LFB 为 月 球 借 力 转移 。 

表 1 EMI2 点 Halo 轨道 不 同 转移 方式 比较 


Table 1 Comparison of several types of transfers to 


the EML2 Halo orbits 


> 


参考 。 Halo 轨道 停泊 轨 动力 学 轨道 转 
文献 幅 值 /km 道 高 度 /km 模型 移 方式 
日 8000 200 CR3BP IT 
图 53000 185 JPL WSB 
图 53000 185 CR3BP WSB 
加 6000 200 CR3BP LFB 
回 5000 200 JPL LFB 
内 5000 400 CR3BP LFB 
8] 5000 200 CR3BP LFB 
9] 8000 167 BCRFBP WSB 


本 文 针对 以 上 问题 ,结合 
性 ,利用 遗传 算法 与 微分 修正 算法 结合 的 混合 优化 
技术 ,设计 了 同时 考虑 地 球 停泊 轨道 高 度 、 倾 角 、 航 


平 动 点 不 变 流 形 的 特 


总 速度 增 量 /(ms"') 


1 总 速度 增 量 ( 直方 图 ) 与 飞行 时 间 ( 曲线 ) 


Fig.1 Total velocity increment( bar) and time of flight ( curve) 


迹 角 等 约束 的 燃料 最 优 地 月 段 转移 轨道 。 同 时 ,在 


不 同 应 用 背景 下 ,通过 微分 修正 算法 调整 稳定 流 形 ， 
得 到 了 满足 月 球 借 力 约束 的 月 球 至 Halo 轨道 的 转 
移 轨道 。 重 点 分 析 了 在 不 同月 球 借 力 约束 下 ,不 同 
幅 值 Halo 轨道 的 不 同 入 轨 点 与 飞行 时 间 、 机 动 速度 
增 量 之 间 的 关系 , 总 结 相 应 的 变化 规律 , 最终 确定 
HOI 点 的 选择 策略 ,对 以 EML2 点 Halo 轨道 为 目标 
轨道 的 深 空 探测 器 转移 轨道 设计 具有 借鉴 意义 。 


1 转移 轨道 模型 及 设计 方案 分 析 


1.1 圆 型 限制 性 三 体 模 型 

针对 地 球 停泊 轨道 与 EML2 点 Halo 轨道 之 间 
的 转移 轨道 ,本 文 基于 圆 型 限制 性 三 体 模型 进行 分 
析 , 则 在 会 合 坐标 系 下 , 探测 器 的 动力 学 模型 可 表示 
为 [qd : 


2 /aU(n)Y 
| ) +h(v) (1) 
式 中 : 
[x -QD (+ Ei 
aU(r) _ yy py 
or eo 
ek 必 
h(y) = By，- 2xz;0] ,7,7 分 别 为 探测 器 与 地 球 


和 月 球 间 的 距离 ,地 - 月 系 中 ,人 ~ 1.215 x 10”， 
kk,y,z,*,] 表征 探测 器 在 会 合 坐标 系 中 的 状态 


[= 


星 。 


1.2 Halo 轨道 及 状态 转移 矩阵 
EMI2 点 是 位 于 月 球 背 面 的 共 线 平 动 点 , 其 附 
近 存 在 的 周期 Halo 轨道 能 够 为 执行 特殊 空间 任务 
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的 探测 器 提供 理想 航行 位 置 。 本 文 针 对 不 同 幅 值 
Halo 轨道 进行 研究 "和 , 即 对 法 向 幅 值 由 2000 km ~ 
10000 km 的 9 组 参数 进行 分 析 。 

状态 转移 矩阵 更 ( t,to) 用 于 描述 下 一 时 刻 目标 
点 的 状态 改变 量 与 当前 时 刻 的 状态 变化 量 之 间 的 关 
系 。 在 圆 型 限制 性 三 体 问题 模型 中 , 状态 转移 矩阵 
满足 以 下 微分 方程 。 

Bi,to) = A(i) B(t,to) (2) 
式 中 : 状态 转移 矩阵 的 初 值 为 6 x6 的 单位 矩阵 , 即 
$B(4o,to) = 和 ;矩阵 4( 为 非 线 性 系统 的 雅 可 比 
和 矩阵 , 满足 
03.3 六 


a 2 0 
U, U, VU, -200 
U. Us: U, 0 0 0 


U0 为 U(r) 关于 r = kk,y 妆 "的 二 阶 偏 导数 , /为 
单位 矩阵 。 联 立 求解 式 ( 1) 与 式 ( 2) ,能 够 确定 任意 
时 刻 的 状态 转移 矩阵 。 
1.3 ”转移 轨道 设计 方案 
本 文通 过 脉冲 修正 EML2 点 Halo 轨道 的 稳定 流 
形 并 结合 月 球 借 力 技术 ,给 出 了 三 脉冲 转移 轨道 设计 
方案 。 同时 ,为 了 实现 对 月 球 背面 近 距 离 观 测 与 未 来 
为 月 球 背面 搭建 无 线 电 望 远 镜 时 进行 物质 运输 投放 
等 任务 ,设计 了 短 时 间 绕 月 飞行 的 四 脉冲 转移 轨道 。 
停泊 轨道 


要 
> 
Eg 


转移 轨道 


Halo 轨 道 


逃逸 点 “AFano 
(a) 三 脉冲 设计 方案 
(a) Design scheme using three-impulse 


停泊 轨道 


由 
f: 


转移 轨道 


逃逸 点 LEO 
(b) 四 脉冲 设计 方案 
(b) Design scheme using four-impulse 
2 ”转移 轨道 设计 方案 


Fig.2 Transfer orbit design scheme 
在 三 脉冲 与 四 脉冲 设计 方案 中 , 重点 分 析 了 不 
同月 球 借 力 约束 条 件 下 ,不 同 幅 值 Halo 轨道 的 不 同 
入 轨 点 对 转移 轨道 参数 的 影响 ,最 终 确 定 不 同 设计 


方案 的 HOI 点 选择 策略 ,如 图 2 所 示 。 
2 约束 条 件 与 分 析 方 法 


为 了 利用 与 Halo 轨道 相连 的 稳定 流 形 ,本 文采 
用 逆向 积分 策略 设计 地 月 平 动 点 转移 轨道 。 在 考虑 
多 约束 条 件 下 ,分 别 对 月 球 至 Halo 轨道 与 地 - 月 间 
的 转移 轨道 进行 分 析 。 
2.1 月 球 至 目标 Halo 轨道 段 

在 转移 轨道 设计 过 程 中 ,当月 球 借 力 点 考虑 不 
同 约束 时 ,将 对 速度 增 量 等 参数 产生 较 大 的 影响 。 
为 了 分 析 不 同 约束 条 件 下 的 借 力 效 果 , 首 先 定义 月 
球 借 力 约束 , 即 近 月 高 度 h, 、 航 迹 角 y,, 、 借 力 方位 
角 6, ,其 中 , 航 迹 角 y, 为 速度 矢量 和 垂直 于 位 置 矢 
量 的 法 线 之 间 的 夹 角 ; 方位 角 6, 为 月 球 借 力 点 轨道 
平面 与 会 合 系 中 xy 平面 所 形成 的 二 面 角 。 航 迹 角 
与 方位 角 的 描述 如 图 3 所 示 。 


3 航 迹 角 y，( 左 ) 与 方位 角 5,，( 右 ) 约束 


Fig.3 Constraints of track angle y,, ( left) and 


azimuth 6,, ( right) 


假设 会 合 系 中 借 力 点 位 置 、 速 度 矢量 为 : rw = 
bi 373zal“ 与 VV = i,, ;76 326] ”; 则 探测 器 相 
对 于 月 球 质心 的 位 置 矢 量 表示 为 Ps, = Ekk;, -1+ 
symszml ”， 月 球 半 径 为 R, ,月 球 借 力 约束 条 件 满 
足 式 (3) , 即 
所 VC 一 工 + AL) + yh, 十 2 一 民 。 


Ym = arcsin[ (xm -1+ 1) Xm 二 mym 浊 | 
|P lv, 
0,, = arccos| (x6 ~ 1 +H) ym — | 
|g. | 
(3) 
inzm pa 


式 中 : Os, 一 本 (xm 人 1 + 1) i 7 文 中 


(xm 1 + 人 Ya Yea 


| 为 2 东 6 数 。 
由 于 目标 Halo 轨道 的 稳定 流 形 能 够 延伸 到 月 
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球 附近 ,可 利用 稳定 流 形 为 转移 轨道 的 设计 提供 初 
值 。 同 时 , 在 给 定 约束 条 件 下 ,微分 修正 作为 一 个 局 
部 优化 算法 ,具有 快速 收敛 特性 。 因 此 ,采用 最 小 二 
乘 微分 修正 算法 对 流 形 进行 脉冲 调整 , 使 得 转移 轨 
道 未 端 满足 期 于 的 月 球 借 力 约束 ,具体 步骤 如 下 : 

Step 1. 针对 转移 轨道 末端 约束 条 件 , 本 文选 定 
控制 变量 Ce 为 HOI 点 三 轴 速 度 矢量 VV 与 飞行 时 间 
7 。 通 过 上 比较 约束 条 件 的 实际 值 与 期 望 值 C,, , 确 
定 约束 矢量 F( C,) , 即 


Co 区， 太太， 如 Bo 
hh Ss ha (4) 
F(C,) = | siny。- sinya, 


COs6,, — cosO。 


moon 


Step 2. 结合 状态 转移 和 矩阵 $B ,推导 约束 矢量 
相对 于 控制 变量 的 微分 关系 式 。 具 体 地 , 根据 链 规 


则 ,约束 矢量 满足 
aF aF aF, oF, 
OH， = yo: + Vo + es + roT 
(0. 韦 6.j5Y 击 (EV, + tar, 37 
ores, 9 Vs, Ors, 9 Vi, 
(5) 
ee cn 
: 一 ,一 一 的 过 程 Ne 
式 中 :9 的 推导 过 程 如 式 (6) 所 示 ， 37 5 
oF RR i a 
与 ap 的 表达 式 类 似 ; 6F, 为 式 (4) 中 约束 矢量 的 


第 ; 行 表 达 式 ; @. ,更 .为 状态 转移 矩阵 中 3 x3 子 


D, 9, > > es 
和 矩阵, 即 态 = | | ; an 为 月 球 借 力 点 加 速度 
0D 二、 
矢量 。 
arm 
01， oF, 9Xs, | | oV, | _ 
OV 6Xm OV ors, avi | ay 
oV, 
[要 |g 
Or, l 9Vi, 9D, 
Om 
oF, oF, 9X a oT | 
07 Xe 67 bors, Ye ov 
oT 
| | 
ar avo, la, 


由 式 (5) 可 知 ,求解 2 


ot 
于 控制 变量 的 微分 关系 式 。 
(1) 借 力 高 度 约束 
结合 式 (3) ` 式 (4) 及 矢量 求 导 法 则 ,月 球 借 力 
点 高 度 约束 的 微分 关系 式 满足 
oF ooh, Pr 


即 可 确定 约束 矢量 关 


or Pl | 
三 0 .3 
fm 


(2) 飞行 航 迹 角 约 束 
探测 器 相对 于 月 球 的 飞行 航 迹 角 约束 的 微分 关 


1 Ve Pi 
or Pa Teo 
Es 
Tall 1 EY 


(3) 借 力 方位 角 8, 约束 
普 力 方位 角 6,, 约束 关于 位 置 、 速 度 矢量 的 微分 


oF, 1 (Ps,V,) VE EE | Vs, Ps, 


> 一 Kk + CoSO ， 
Ors, | Or, | : | Os, 下 
oF (Po Pe lB DY 
= 1 kK, + 站 | | cos8 ， 
aVi, oul lg 


(9) 
式 中 : 
Kk = bs, 一 xm30] ,k= EF- yra? Xen 一 + 人 ;0] o 
Step 3. 在 微分 修正 过 程 中 , 选 定 积分 终止 条 件 
为 零 航 迹 角 ,可 以 保证 转移 轨道 较 好 地 满足 借 力 点 
约束 条 件 。 因 此 ,处 理 约束 条 件 时 , 设 定 f= siny， 
Ss Sinydm = 0 ,由 式 (5) 可 推导 67 = 67( 6V) ;代入 
普 力 高 度 与 方位 角 约 束 的 微分 关系 式 中 ,满足 
oF, 9F./9T oF, 
[551] _ | am ~ oF,/9T ov a 
OF, oFs 9F;/90T oF, 
| aV' oF,/9T ov 
利用 最 小 二 乘法 迭代 求解 式 ( 10) ,可 确定 6V ， 
进而 修正 HOI 点 三 轴 速 度 矢量 ,使 得 最 终 的 转移 轨 
道 末 端 收敛 于 期 望 的 约束 值 。 
2.2 ”地 球 至 月 球 转移 轨道 段 
为 了 体现 圆 形 地 球 停泊 轨道 特性 ,考虑 末端 约 
束 条 件 包括 轨道 高 度 h,、 轨 道 倾 和 角 i, 、 升 交点 赤 经 
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人 .及 地 - 月 转移 飞行 时 间 7, ,并 以 相对 地 球 的 航 迹 
角 y, 为 积分 终止 条 件 。 

针对 地 - 月 转移 轨道 特性 研究 ,本 文采 用 遗传 
算法 ( Genetic algorithm，GA) 与 微分 修正 算法 
( Differential correction，DC) 相 结 合 的 混合 优化 方 
法 。 一 方面 ,充分 利用 遗传 算法 的 全 局 收敛 特 性 与 
微分 修正 的 局 部 快速 收敛 特 性 ; 另 一 方面 ,遗传 算法 
的 优化 结果 可 作为 微分 修正 算法 的 初 值 ,避免 了 微 
分 修正 的 初始 猜想 。 基 于 混合 优化 算法 , 设计 了 燃 
料 最 优 的 地 月 转移 轨道 。 

(1) 遗传 算法 全 局 优化 

遗传 算法 是 基于 基因 遗传 学 原理 与 自然 选择 策 
略 的 智能 优化 算法 " 。 利 用 群体 搜索 策略 与 个 体 
之 间 的 信息 交换 ,以 适应 度 值 判断 个 体 的 好 坏 , 并 结 
合 选择 、 交 叉 及 变异 等 一 系列 遗传 操作 , 最 终 运 近 问 
题 的 最 优 解 。 

在 地 月 转移 轨道 设计 时 ,优化 变量 为 借 力 点 速 
度 增 量 AV,。= ”人 V,,AV, ;AV ,以 燃料 最 优 为 性 
能 指标 ,分 析 可 知 , Halo 轨道 入 轨 速 度 增 量 对 总 速 


度 增 量 的 影响 较 小 。 因 此 ,选取 目标 函数 J 为 地 
球 逃 逸 点 与 月 球 借 力 机 动 点 速度 增 量 大 小 , 满足 
Jea = AT + AVoo, (11) 
式 中 : AV,,,, 为 月 球 借 力 点 三 轴 速 度 增 量 大 小 , 即 
AV, 上 。 优 化 过 程 中 同时 考虑 地 球 停泊 轨道 高 度 、 
轨道 倾角 、 升 交点 赤 经 及 飞行 时 间 约 束 , 定 义 为 不 等 
式 约束 ,满足 
hi < h, < ho, 
(12) 


(2 < 1, < OQ 
Ti < 了 < Th 
式 中 : 下 标 min 与 max 分 别 表 示 在 遗传 算法 优化 中 ， 
设 定 的 约束 条 件 最 小 值 与 最 大 值 。 
( 2) 微分 修正 局 部 优化 
通过 遗传 算法 能 够 使 地 月 转移 轨道 末端 位 于 期 
望 的 地 球 停泊 轨道 附近 ,并 且 轨 道 根 数 满足 我 国 发 
射 要 求 。 因 此 ,在 微分 修正 局 部 优化 中 ,只 考虑 轨道 
高 度 与 航 迹 角 约 束 。 


h, -hh 
we = | ee | (13) 
siny, — Siny ag 
同 理 , 结 合式 (7) .(8) 与 式 ( 10) 的 微分 关系 ,能 


够 设计 出 同时 满足 轨道 高 度 与 航 迹 角 约 束 的 地 - 月 
转移 轨道 。 


综 上 所 述 ,混合 优化 求解 过 程 如 图 4 所 示 。 


定义 轨道 设计 参数 AV =[AV,AV ,AV], 
给 定 设计 参数 上 下 限 区 问 


确定 月 球 供 力 点 的 轨道 初 什 
已 [yz X+A 及 HA Vz+ A 了 


产生 新 的 轨 
道 设计 参数 


由 设计 参数 确定 轨道 
初 值 ， 并 结合 动力 学 


达 代 修 
六 程 进行 积分 过 代 修正 


控制 变量 
AD, 


多 轨道 末端 
满足 约束 条 件 


4 混 口 合 优化 算法 流程 图 
Fig.4 Flow chart of the hybrid optimization algorithm 


3 仿真 结果 分 析 


本 文 以 地 月 系统 L2 点 附近 不 同 法 向 幅 值 的 
Halo 轨道 为 目标 轨道 ,在 考虑 不 同月 球 借 力 约束 与 
地 球 停泊 轨道 约束 条 件 下 , 通过 脉冲 调整 稳定 流 形 
并 结合 混合 优化 算法 ,设计 了 三 脉冲 与 四 脉冲 的 平 
动 点 转移 轨道 。 

为 充分 分 析 目 标 Halo 轨道 上 不 同 HOI 点 对 飞 
行 时 间 与 机 动 速度 增 量 的 影响 ,首先 将 Halo it 
照 时 间 等 分 为 360 个 点 ,定义 Halo 轨道 与 会 合 系 xy 
平面 相交 的 两 点 中 距离 月 球 较 远 的 点 为 1 点 ,积分 
一 周 后 对 Halo 轨道 上 各 点 进行 编号 。 

3.1 三 脉冲 转移 轨道 方案 分 析 

在 三 脉冲 转移 轨道 设计 过 程 中 , 当 保持 其 他 条 
件 不 变 的 情况 下 , 借 力 高 度 越 低 , 借 力 飞行 的 效果 越 
显著 ”。 因 此 ,本 文 研究 借 力 约束 条 件 对 燃 耗 等 参 
数 9 影响 时 ， 对 固定 的 借 力 高 度 进行 分 析 , 任务 约束 
如 表 2 所 示 。 
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表 2 任务 约束 条 件 


Table 2 The mission constraints 


参数 约束 条 件 
地 球 停 泊 轨 道 高 度 h,/km 200 
偏心 率 。 0 
倾角 范围 (i ;is) /( 9?) (20,60) 
RAAN 范围 (0 ,02%) A(°) ( -70;0) 
地 -月 段 飞行 时 间 ( Tu ， 7) /天 (4,6) 
发 射 时 间 ( 年 /月 /日 ) 2020/01/01 
航 迹 角 y,,y,, A(“) 0 
月 球 借 力 点 高 度 h,, /km 100 
位 置 误差 /m 1 
角度 误差 /(“) 1x10-3 


3.1.1 只 考虑 借 力 高 度 约束 的 借 力 效果 分 析 

在 表 2 的 任务 约束 下 , 仅 施 加 月 球 借 力 高 度 约 
束 ,分 析 不 同 法 向 幅 值 Halo 轨道 .不 同 HOI 点 与 各 
个 机 动 点 速度 增 量 大 小 及 飞行 时 间 的 关系 。 由 
5 ~7 可 知 : ( 1) 只 考虑 借 力 高 度 约束 时 , 随 着 Halo 
轨道 法 向 幅 值 不 断 增 大 ,HOI 点 空间 位 置 逐步 改变 ， 
为 了 满足 借 力 高 度 约束 ,需要 消耗 更 大 的 速度 增 量 
用 于 调整 稳定 流 形 。 因 此 ,针对 相同 编号 的 HOI 点 ， 
不 同 幅 值 Halo 轨道 所 需 总 速度 增 量 呈现 增长 趋势 且 
变化 显著 ; (2) 针对 各 个 机 动 点 施加 的 速度 增 量 大 小 ， 
Halo 轨道 入 轨 速 度 增 量 单调 变化 , 法 向 幅 值 对 HOI 
点 速度 增 量 的 影响 较 小 。 同 时 ,月 球 借 力 点 与 地 球 逃 
逸 点 速度 增 量 是 影响 总 速度 增 量 差异 的 主要 因素 ; 
(3) 相 比 于 地 月 段 飞 行 时 间 , 月球 至 Halo 轨道 段 耗 时 
较 长 ,与 总 飞行 时 间 变 化 一 致 ,保持 单调 递增 。 


3500 


3480 


N46000 km | 
4.=5000 km  ] 


4_=4000 km 


RL km 


30 60 90 120 150 180 210 240 270 300 330 360 
Halo 轨 道 入 轨 点 


5 不 同 幅 值 Halo 轨道 HOI 点 消耗 的 总 速度 增 量 
Fig.3 Total speed increment at HOI of Halo orbits 


with different amplitude 


结合 图 5 与 图 6, 随 着 Halo 轨道 法 向 幅 值 逐步 
增 大 ,仅仅 考虑 月 球 借 力 高 度 约束 并 不 能 有 效 地 降 
低 借 力 点 的 速度 增 量 。 通 过 分 析 借 力 方位 角 6, 的 
变化 ,如 图 8 所 示 , 当 Halo 轨道 法 向 幅 值 为 2000 km 
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Fig.7 Time of flight in each transfer orbit segment 


时 , 借 力 方位 角 6,,，e (4°,5°) , 随 着 法 向 幅 值 的 堆 
大 ,相同 HOI 点 对 应 的 借 力 方 位 角逐 渐 增 加 。 
此 ,本 文 在 借 力 高 度 与 航 迹 角 约 束 的 基础 上 ,进一步 
添加 方位 角 约 束 ,分 析 不 同方 位 角 对 速度 增 量 与 总 
飞行 时 间 的 影响 ,进而 确定 降低 飞行 燃 耗 的 策略 。 


SS 
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Fig.8 Gravity-assist azimuth for Halo orbits 


with different amplitude 
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3.1.2 借 力 方位 角 对 借 力 效果 的 影响 分 析 

通过 脉冲 调整 稳定 流 形 可 确保 轨道 末端 满足 借 
力 高 度 与 航 迹 角 约 束 ,然而 , 当 添 加 方位 角 约 束 后 , 音 
分 稳定 流 形 不 易 收敛 到 期 望 的 借 力 点 。 因 此 ,以 下 针 
对 不 同方 位 角 约 束 ,分 析 不 同 幅 值 Halo 轨道 的 转移 


轨道 特性 ,数据 对 比 后 确定 HOI 点 的 选择 策略 。 

9 描述 了 不 同 借 力 方位 角 约 束 下 飞行 时 间 - 
速度 增 量 曲线 , 左 图 中 线段 间 空 白 区 域 表征 稳定 流 
形 未 能 较 好 地 满足 约束 条 件 , 右 图 为 总 飞行 时 间 
7 e (18,24) 天 的 放大 图 。 


HOI 点 


仅 含 高 度 约束 
250 r ' 一 5- 
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9 4. =5000 km 时 不 同 借 力 方位 角速度 增 量 曲线 


Fig.9 Speed increment curves for different gravity-assist azimuth when 4. =5000 km 


(1) 针对 相同 HOI 点 , 随 着 借 力 方位 角 的 减 小 ， 
佛 力 飞行 对 入 轨 点 速度 增 量 的 影响 逐渐 增强 ,但 借 
力 点 所 需 的 速度 增 量 呈现 下 降 趋 势 , 总 速度 增 量 相 
比 于 仅 考 虑 借 力 高 度 约束 时 不 断 降低 ,最 大 差 值 为 
58.52 m/s。 

(2) 针对 不 同 HOI 点 ,选取 借 力 方位 角 5,，= 4° 
分 析 , 在 7e (12,17) U (18,24) 天 区 间 内 ,总 速度 
增 量 存在 两 个 极 小 值 , 即 Node = 143, 耗 时 14. 81 
天 ,总 燃 耗 3329. 94 m/s; Node =310, 飞行 时 间 与 总 
速度 增 量 分 别 为 21. 77 天 ,3322. 65 m/s。 同 时 , 满 
足 约 束 的 入 轨 点 数目 随 着 借 力 方位 角 的 变化 而 减 
综合 分 析 , 较 好 地 满足 借 力 约束 的 HOI 点 选取 
范围 为 B0,160) 与 B40,315) 。 

表 3 选取 幅 值 为 5000 km 的 Halo 轨道 ,对 比分 
析 不 同 借 力 约束 下 ,机动 点 消耗 的 速度 增 量 大 小 ,其 
中 各 个 HOI 点 的 第 一 行 数据 为 仅 考 虑 借 力 高 度 与 
航 迹 角 约 束 的 情况 。 

10 与 图 11 表明 : (1) 当 Halo 轨道 法 向 幅 值 
增 大 时 , 总 速度 增 量 逐步 增加 , 满足 Ar e 
(3320 m/s,3400 m/s) 。 进 一 步 说 明了 大 幅 值 Halo 


轨道 的 稳定 流 形 在 月 球 处 的 借 力 效果 较 差 , 必须 施 
加 更 大 的 脉冲 对 流 形 进 行 修正 以 满足 给 定 的 约束 ; 
(2) 在 给 定 的 两 个 入 轨 点 区 间 内 , 随 着 入 轨 点 的 变 
化 ,HOI 点 速度 增 量 总 体 呈 现 先 减 小 后 增 大 的 趋势 ， 
而 借 力 点 增 量 变化 相反 。 
表 3 不 同 借 力 约束 条 件 下 参数 比较 
Table 3 Parameter comparison under different 


gravity-assist constraint conditions 


借 力 广 HOI 点 机 动 借 力 点 机 动 三 脉冲 总 
HOI 点 位 角 /(”) 速度 境 量 速度 增 量 速度 增 量 
/(m*s-!) /(m*s-!) /(m*s-!) 
0.51 18.34 226. 16 3372. 88 
90 4.00 38.64 189.39 3348.82 
8.00 20.71 203.15 3345.11 
0.97 9.33 233.33 3372.68 
120 4.00 16.92 193. 82 3331.85 
8.00 11.04 206. 32 3339.52 
2.05 0.22 248.39 3381.32 
270 4.00 13.15 190. 85 3324.48 
8.00 6.50 202.73 3330.77 
2.04 0.03 248.38 3381.11 
300 4.00 12.05 190. 52 3323. 00 
8.00 5.28 204.01 3331.04 
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Fig. 11 


表 4 给 出 了 不 同 幅 值 条 件 下 对 应 的 最 佳 入 轨 
点 、 停 泊 轨 道 等 各 项 参数 。 

综 上 所 述 ,选取 幅 值 4. = 10000 km 的 Halo 轨 
道 为 目标 轨道 ,设计 了 三 脉冲 转移 轨道 ,如 图 12 所 
示 , 其 中 虚线 为 只 考虑 借 力 高 度 与 航 迹 角 约 束 , 实 线 
为 添加 方位 角 约 束 的 转移 轨道 。 在 借 力 高 度 与 航 迹 
角 约 束 条 件 下 ,总 飞行 时 间 为 13. 25 天 , 地 球 逃 逸 


Speed increment at HOI and gravity-assist points when 6 


=6° 


m 


点 、 月 球 借 力 点 与 HOI 点 速度 增 量 大 小 分 别 为 
3133.41 m/s,685. 01 m/s 与 9. 67 m/s。 当 添加 借 
力 方位 角 约 束 后 ,三 处 机 动 点 速度 增 量 大 小 分 别 为 
3120. 33 m/s,190. 24 m/s 与 39. 62 my/s, 总 飞行 时 
间 变 化 较 小 , 即 13. 63 天 。 比 较 两 种 情况 可 知 , 在 飞 
行 时 间 相 近 条 件 下 ,考虑 借 力 方 位 角 后 的 总 速度 增 
量 显 著 降低 ,相差 477. 9 m/s。 
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表 4 方位 角 6,，= 4° 时 不 同 幅 值 的 借 力 效果 对 比 
Table 4 Comparison of gravity-assist effects in 


different amplitude when 6, = 4° 


m 


Halo 轨道 ” 最 佳 入 倾角 / 升 交 ”飞行 时 间 ”总 速度 堆 
幅 值 /km ” 轨 点 HOI ， 点 赤 经 /( ") /天 ” 量 /(m*s-1!) 
S5000 143 35.8/—10.9 14.81 3329.94 
6000 141 23.0/ -19.6 14.89 3333.28 
7000 135 23217 三 755 14.57 3337.29 
8000 128 21.1/ -21.1 14.08 3341. 80 
9000 126 21.1/ -26.5 13.96 3345.74 
10000 120 21.6/ -35.3 13.63 3350.19 

0.6 
合 0.4 
02 加 
所 0 佛 力 高 度 与 航 迹 角 约束 
Ka 


1 
x/(105 km) 


图 12 不同 借 力 约束 条 件 的 三 脉冲 转移 轨道 


Fig.12 Threeimpulse transfer trajectories under different 


gravity-assist constraint conditions 


3.2 ”四 脉冲 绕 月 转移 轨道 方案 分 析 

由 于 月 球 背 面 在 未 来 的 空间 研究 中 具有 巨大 的 
应 用 价值 ,如 何在 燃 耗 较 低 的 条 件 下 ,实现 探测 器 对 
月 球 背 面 近 距 离 观 测 。 针 对 此 问题 ,采用 与 三 脉 ; 
设计 方案 类 似 的 方法 , 设计 了 短 时 间 绕 月 飞行 的 四 
脉冲 转移 轨道 ,分 析 了 不 同 HOI 点 在 不 同 借 力 方位 
角 约 束 下 的 转移 轨道 特性 。 目 标 函 数 为 

Jes =APF+AY + AV,, 
式 中 : AV,, , AV, 分 别 表示 月 球 两 次 机 动 点 的 速度 
增 量 大 小 。 

13 中 ,探测 器 以 3118.35 m/s 的 速度 增 量 从 
轨道 高 度 为 200 km、 倾 角 24°、 升 交点 赤 经 - 31 的 
地 球 停泊 轨道 逃逸 ,飞行 4. 42 天 后 到 达 高 度 为 
400 km 的 月 球 借 力 机 动 点 1, 施加 脉冲 250. 92 m/s 
形成 1. 33 天 的 绕 月 飞行 轨道 , 当 飞 行 高 度 为 
100 km ,方位 角 为 4? 时 进行 第 二 次 借 力 机 动 进入 月 
球 - Halo 轨道 段 , 增 量 大 小 为 51. 38 m/s, 航行 
8.97 天 后 ,以 15.61 m/s 的 增 量 完成 入 轨 操 作 。 

14 与 表 5 进一步 分 析 了 四 脉冲 轨道 设计 方 


案 中 ,不 同 借 力 约束 与 机 动 点 速度 增 量 之 间 的 关系 ， 
表 5 中 各 项 法 向 幅 值 包含 的 四 组 数据 依次 对 应 方位 
角 6,, = 4°,6°,8°,10* 。 分 析 可 知 : 随 着 借 力 方位 角 
的 减 小 ,满足 约束 的 HOI 点 数目 及 速度 增 量 、 借 力 
点 1 速度 增 量 、 总 速度 增 量 的 变化 趋势 与 三 脉冲 设 
计 方 案 结果 相似 ,而 借 力 效果 对 飞行 时 间 、 借 力 点 2 
9 影响 较 小 。 


速度 增 量 


z/(10km) 


一 1 
nOIOD 


1 
-5 0 x(10km) 
13 4. =5000 km 的 四 脉冲 转移 轨道 
Fig. 13 Fourimpulse transfer trajectory when 4. =5000 km 
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100 130 140 160 180 18 280 30 360 


Halo 轨 道 入 轨 点 Halo 轨 道 入 轨 点 
14 4. =5000 km 时 不 同 借 力 方位 角 约 束 分 析 


Fig.14 Analysis of constraints of different gravity-assist 


azimuth when 4. =5000 km 


表 5 Node =125 时 不 同方 位 角 转 移 轨道 参数 对 比 
Table 5 Parameter comparison of transfer trajectory with 


different azimuth when Node = 125 


Halo 轨 借 力 点 1 借 力 点 2 HOI 点 机 ”四 脉冲 总 
道 法 向 速度 增 量 。 速度 增 量 ” 动 速度 增 量 ”速度 增 量 
幅 值 /km /ms-) /ms-) /ns /(m*s-!) 
251.29 51.38 16.24 3437.33 
256. 86 51.38 12. 92 3440.76 
5000 
265.67 51.37 10.30 3448. 80 
281.40 51.34 8.73 3466.27 
249.77 51.40 40.81 3460. 16 
258.56 51.46 34.12 3464.17 
10000 
270.35 51.48 29.74 3474.09 
292.76 51.48 25:77 3497. 35 
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4 结 论 


本 文 针对 地 月 L2 点 不 同 任务 需求 的 转移 轨道 
设计 问题 ,基于 不 变 流 形 理论 与 混合 优化 技术 , 重点 
研究 了 不 同月 球 借 力 约束 与 HOI 点 对 转移 轨道 飞 
行 时 间 、 机 动 速 度 增 量 的 影响 。 研 究 结果 表明 ,在 圆 
型 限制 性 三 体系 统 中 : 

(1) 针对 小 幅 值 Halo 轨道 , 仅 考虑 借 力 高 度 与 
航 迹 角 约 束 即 可 设计 出 满足 要 求 且 燃 耗 较 低 的 三 肪 
冲 转移 轨道 。 其 中 总 速度 增 量 主要 取决 于 借 力 点 与 
地 球 逃 逸 点 所 施加 的 脉冲 大 小 ,总 飞行 时 间 随 着 
HOI 点 的 增 大 呈现 单调 变化 。 

(2) 针对 固定 幅 值 Halo 轨道 分 析 , 添加 借 力 方 
位 角 约 束 后 , 满足 要 求 的 HOI 点 数目 逐渐 减少 ,但 
相 比 于 仅 考 虑 借 力 高 度 与 航 迹 角 约 束 , 月 球 借 力 机 
动 增 量 与 总 速度 增 量 显著 下 降 , 而 借 力 方位 角 的 变 
化 对 飞行 时 间 的 影响 较 小 。 

(3) 针对 不 同 幅 值 Halo 轨道 分 析 , 选 定 借 力 方 
位 角 和 与 Halo 轨道 入 轨 点 后 ,由 于 流 形 特性 的 影响 ， 
随 着 法 向 幅 值 的 增 大 , HOI 点 施加 脉冲 依次 增加 , 借 
力 点 速度 增 量 降低 ,总 燃 耗 不 断 增 大 。 同 时 ,在 四 脉 
冲 设计 方案 中 , 借 力 点 2 所 需 速 度 增 量 变化 较 小 。 

以 上 的 研究 与 分 析 , 对 于 考虑 月 球 借 力 的 地 月 
平 动 点 转移 轨道 的 设计 与 应 用 、 参 数 选 择 等 具有 重 
要 的 借鉴 意义 。 
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